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ABSTRACT 

We construct relativistic equilibrium models of differentially rotating neutron stars and show that 
they can support significantly more mass than their nonrotating or uniformly rotating counterparts. 
We dynamically evolve such "hypermassive" models in full general relativity and show that there do 
exist configurations which are dynamically stable against radial collapse and bar formation. Our results 
suggest that the remnant of binary neutron star coalescence may be temporarily stabilized by differential 
rotation, leading to delayed collapse and a delayed gravitational wave burst. 
Subject headings: black hole physics — relativity — stars: neutron — stars: rotation 



1. INTRODUCTION 

One of the most important characteristics of a neutron 
star is its maximum allowed mass. The maximum mass is 
crucial for distinguishing between neutron stars and black 
holes in compact binaries and in determining the outcome 
of many astrophysical processes, including supernova col- 
lapse and the merger of binary neutron stars. 

Observations of binary pulsars suggest that the individ- 
ual stars in such systems have masses very close to 1.4 
Mq (Thorsett et al. 1993). If mass- loss during the final 
binary coalescence can be neglected, the remnant of such a 
merger will then have a rest mass exceeding 3 Mq . If this 
mass is larger than the maximum allowed mass for neu- 
tron stars, then the merger will lead to prompt collapse 
to a black hole on a dynamical timescale (ms). If, how- 
ever, neutron stars can support such a high mass, at least 
temporarily, then the merger may result in a high mass, 
quasi-equilibrium neutron star, which only later may col- 
lapse to a black hole. The two different outcomes may have 
important consequences for gravitational wave signals and, 
possibly, gamma ray burst models. 

The maximum mass of a cold, nonrotating, spheri- 
cal neutron star is uniquely determined by the Tolman- 
Oppenheimer-Volkoff equations and depends only on the 
cold equation of state. For most recent equations of state 
this maximum mass is in the range of 1.8 - 2.3 Mq (Akmal, 
Pandharipande and Ravenhall 1998), significantly smaller 
than the mass expected for the remnant of a binary neu- 
tron star merger. 

Thermal pressure and uniform rotation can provide ad- 
ditional support, and may stabilize slightly more massive 
stars. In this paper, we point out that differential rota- 
tion can significantly increase the maximum allowed mass 
of neutron stars and may temporarily stabilize the rem- 
nant of binary neutron star mergers. Recent fully relativis- 
tic simulations of binary neutron star mergers (Shibata & 
Uryu 1999) show that merger remnants are indeed differ- 
entially rotating (as suggested by several Newtonian sim- 
ulations) and that they may support masses much larger 
than the maximum allowed mass of spherical stars. 



In the case of a head-on collision from infinite separa- 
tion, thermal pressure alone may support the merged rem- 
nant of progenitors (Shapiro 1998). Thermal pressure is 
likely to have a much smaller effect for coalescence from 
the innermost stable circular orbit since shock heating on 
impact is less pronounced, and will be dissipated by neu- 
trino emission in ^ 10 s. 

Rotation can further increase the maximum allowed 
mass. The maximum mass of a uniformly rotating star is 
determined by the spin rate at which the fluid at the equa- 
tor moves on a geodesic and any further speedup would 
lead to mass shedding. This maximum mass can be deter- 
mined numerically and is found to be at most ^ 20% larger 
than the nonrotating value (e.g. Cook, Shapiro, & Teukol- 
sky, 1992, hereafter CST; 1994, and references therein). It 
is therefore unlikely that uniform rotation could support 
the remnant of a binary neutron star merger. Rotating 
equilibrium configurations with rest masses exceeding the 
maximum rest mass of nonrotating stars constructed with 
the same equation of state are referred to as "supramas- 
sive" stars (CST). 

The merger of a binary neutron star system, however, 
will not result in a uniformly rotating object, especially 
since the neutron stars are likely to be close to being irro- 
tational before merger (Bildsten & Cutler 1992; Kochanek 
1992). The remnant is likely to be differentially rotat- 
ing (see Rasio & Shapiro 1999 for discussion and refer- 
ences; Shibata & Uryu 1999). The star's core may then 
rotate faster than the envelope, and it is easy to imagine 
that such a star could support a significantly larger mass 
than its uniformly rotating counterpart (see also Ostriker, 
Bodenhcimer & Lynden-Bell, 1966, where this effect was 
demonstrated for white dwarfs). We refer to equilibrium 
configurations with rest masses exceeding the maximum 
rest mass of a uniformly rotating star as "hypermassive" 
stars. 

In contrast to the maximum mass of the nonrotating 
and uniformly rotating stars, the maximum mass of dif- 
ferentially rotating stars cannot be uniquely defined, since 
the value will depend on the chosen differential rotation 
law. In principle, one might even construct an extensive 
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Keplerian disk around the equator of the star, possibly in- 
creasing the mass of the star by large amounts. Instead 
of constructing such extreme configurations, we seek to 
determine whether a reasonable degree of differential ro- 
tation can have a significant effect on the maximum mass 
of neutron stars. 

Here we adopt a polytropic equation of state and a sim- 
ple rotation law to explore the effects of differential ro- 
tation on the maximum mass. We construct relativistic 
equilibrium models and find that even for modest degrees 
of differential rotation the maximum mass increases sig- 
nificantly, easily surpassing the likely remnant mass of a 
binary neutron star merger. We then evolve high-mass 
models dynamically in full general relativity, and find that 
there do exist models which are dynamically stable against 
both radial collapse and bar formation. These are plausi- 
ble candidates for binary neutron star remnants. 

2. EQUILIBRIUM 

We adopt a polytropic equation of state P = Kp^ , 
where P is the pressure and po the rest-mass density. We 
take the polytropic constant K to be unity without loss of 
generality, and choose the polytropic index n — 1 ^. 

Relativistic equilibrium models of rotating stars have 
been constructed by several authors, including Butter- 
worth & Ipser (1975), Friedman, Ipser & Parker (1986), 
Komatsu, Eriguchi, & Hachisu (1989), CST, Bonazzola et 
al. (1993), and Stergioulas & Friedman (1995). A com- 
parison between several different methods can be found in 
Nozawa et al. (1998). We use the numerical code devel- 
oped by CST, which is based on the formalism of Komatsu, 
Eriguchi, & Hachisu (1989). 

Constructing differentially rotating neutron star models 
requires choosing a rotation law F{Q) = u^u^, where m* 
and u^ are components of the four velocity m", and il is 
the angular velocity. For simplicity we follow CST and 
consider the rotation law F(51) = A'^{Vtc — ri), where Q,c 
denotes the central angular velocity and where the param- 
eter A has units of length. Expressing u* and u^ in terms 
of Vt and metric potentials yields eq. (42) in CST, or, in 
the Newtonian hmit, Vl = fic/(f + Ar'^P sin^ 6). Here we 
have rescaled A and r in terms of the equatorial radius Re'. 
A = A/ Re and f = r/Re- The parameter ^4 is a measure 
of the degree of differential rotation and determines the 
length scale over which fl changes. Since uniform rota- 
tion is recovered in the limit A — > oo, it is convenient to 
parametrize sequences by A^^. 

We construct axisymmetric differentially rotating mod- 
els using a modified version of the scheme adopted in 
CST. Instead of fixing the central density in the iteration 
scheme for each model, we fix the maximum density. This 
change allows us to construct higher mass models in some 
cases, since the central density does not always coincide 
with the maximum density, and hence may not specify a 
model uniquely. Given a value of A, we construct a se- 
quence of models for each value of the maximum density 
by starting with a static, spherically symmetric star and 
then decreasing the ratio of the polar to equatorial radius, 
Rpe — Rp/Re, in small increments. This sequence ends 

^Sincc ftT"/^ has units of length, aU solutions scale aceording to M = 
are physical quantities, and the unbarred quantities are our dimensionless 
*See CST for relativistic definitions of these quantities. 



when we reach mass shedding (for large values of A), or 
when the code fails to converge (indicating the termina- 
tion of equilibrium solutions), or when Rpe — (beyond 
which the star would become a toroid). 
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Fig. 1. — Maximum rest-mass configurations versus maximum 
mass-energy density for differentially rotating n = 1 sequences spec- 
ified by A~^ . Values of /3 = T/\W\ for the models are indicated. The 
mass— density relation for static equilibrium stars (TOV) is shown 
for comparison. Masses in solar masses are calculated by assum- 
ing that the maximum mass for nonrotating stars is 2 Mq. Note 
that even modest differential rotating may easily support > 3Mq, 
the expected mass of binary neutron star merger remnants. We dy- 
namically evolve the model marked with a dot, and show that it is 
dynamically stable (see Fig. 2). 

In Fig. 1 we show the maximum-mass values in each 
sequence as a function of the maximum value of the mass- 
energy density e for different values of A. Even for modest 
differential rotation, we can construct models with masses 
much higher than the maximum mass for static and uni- 
formly rotating stars. Some of these models exceed the 
Kerr limit J /M"^ > 1, where J is the angular momentum 
and M the total mass energy of the star. 

3. STABILITY 

Nonrotating spherical stars are dynamically stable (un- 
stable) against radial modes if dM/dec > {dM/dec < 0) , 
where Cc is the central energy density. The same criterion 
can be applied to sequences of uniformly rotating stars of 
constant J to determine secular stability (Friedman, Ipser, 
& Sorkin 1988). Exact criteria do not exist for the dynam- 
ical stability of rotating stars; however, numerical simula- 
tions of uniformly rotating models suggest that the onset 
of dynamical stability is very close to the onset of secular 
instability (Shibata, Baumgarte, & Shapiro 1999a). 

As an indication of the stability of our models against 
nonaxisymmetric bar-mode formation, we have indicated 
values of the ratio of their kinetic energy T to potential 
energy W, (3 = T/\W\, in Fig. 1 ^. Newtonian stars de- 
velop bars on a dynamical timescale when P ^ /?dyn = 
0.27, while they develop bars on a secular timescale for 

K"/^M, J = K"J, n = K-"/^n, etc, where the barred quantities 
quantities corresponding to K = 1 (compare CST). 



P ^ Ace — 0.14 via gravitational radiation or viscosity 
(Chandrasekhar 1969, 1970; Houser, Centrella, & Smith 
1994). For relativistic stars, (3scc for gravitational wave- 
driven bars is somewhat smaller than for Newtonian stars 
(Stergioulas & Friedman 1998), while /3sec for viscosity- 
driven bars is shghtly larger (Bonazzola, Frieben, & Gour- 
goulhan 1996; Shapiro & Zane 1998). 

t = O.OOP 
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Fig. 2. — Snapshots of contours of the density po""(^s)^'^ of the 
hypermassive star marked with a dot in Fig. 1. We show contours 
of the initial data and after about 3 central orbital periods, both in 
the equatorial plane (left, including the velocity field M'/n*) and in 
a plane containing the z-axis of rotation (right). The star has a rest 
mass of Mo = 0.29, about 60% larger than the maximum nonrotat- 
ing rest mass. The simulation shows that this model is dynamically 
stable against quasi-radial collapse and bar formation. 

To investigate the dynamical stability ol our equilib- 
rium models, we insert them as initial data in a dynamical 
simulation and evolve them in time. We employ a fully 
relativistic code that solves Einstein's equations coupled 
to hydrodynamics in three spatial dimensions plus time 
(Shibata 1999). 

As a candidate for a dynamically stable star, we evolve 
the model with A^^ = 1.0, e^ax = 0.073 and Rpe ~ 0.3 
(marked with a dot in Fig. 1). This model has a rest 
mass about 60% higher than the maximum nonrotating 
rest mass, /3 -- 0.23, Re/M - 5 and J/M'^ ^ 1, and is plot- 
ted in Fig. 2. The orbital period at the equator is about 
three times the orbital period at the center. We show 
contours at t == and after 3.15 orbital periods at the cen- 
ter. Clearly, this model is dynamically stable against both 
quasi-radial collapse to a black hole and bar formation, 
even when small perturbations are included initially. This 
demonstrates that differentially rotating stars can stably 
support significantly higher masses than uniformly rotat- 
ing stars for longer than a dynamical timescale. A more 
systematic study of the dynamical stability of differentially 
rotating neutron stars will be presented in a forthcoming 
paper (Shibata, Baumgarte, & Shapiro 1999b). 

Dynamically stable differentially rotating neutron stars 



are subject to various secular instabilities. The timescale 
for gravitational-wave driven bar-mode formation can be 
estimated from 
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(Friedman & Schutz 1975, 1978), where the average radius 
R and mass M are scaled to values appropriate for a binary 
merger remnant. For /3 ^ 0.2, this yields timescales of 10 
s. The final fate of bar-unstable stars is not known, except 
for incompressible Newtonian spheroids, where in the pres- 
ence of gravitational radiation and viscosity they evolve 
to Jacobian or Dedekind ellipsoids (Chandrasekhar 1969; 
Miller 1974; Shapiro & Teukolsky 1983; Lai & Shapiro, 
1995). Gravitational waves may also drive an r-mode in- 
stability for arbitrarily small rotation rates (see, e.g., Lind- 
blom, Owen & Morsink 1998). For the hot remnants of 
binary neutron star mergers, however, these modes may 
be suppressed by bulk viscosity. 

Magnetic braking and viscosity will eventually bring the 
star into uniform rotation. When a hypermassive star is 
driven to uniform rotation by viscosity or magnetic fields, 
it will undergo catastrophic collapse and/or mass loss. The 
lifetime of a hypermassive star is therefore set by these dis- 
sipative processes. If J /M"^ > 1, angular momentum must 
be dissipated either by radiation or mass loss before the 
star can form a Kerr black hole (cf. Baumgarte & Shapiro 
1998), which may produce a massive, hot and thick disk 
around the newly formed black hole. 

A frozen-in magnetic field will be wound up by differen- 
tial rotation, which may create very strong toroidal fields. 
This process will generate Alfven waves, which can re- 
distribute and even carry off angular momentum. The 
timescale tb for this magnetic braking mechanism is re- 
lated to the Alfven speed va ~ i3/(47rp)^/^ according to 
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Here B is the initial poloidal field along the gradient of 
rj. Strong poloidal magnetic fields can increase the maxi- 
mum allowed mass of neutron stars (Bocquet et al. 1995) 
and contribute to the dissipation of angular momentum 
by dipole radiation, but are subject to a variety of MHD 
instabilities (e.g. Spruit 1999a, 1999b). 

Since the fluid flow in differentially rotating equilibrium 
stars is divergence-free, the viscous timescale tv is deter- 
mined by shear viscosity 
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where r? = Mlp^/^T-"^ (cgs) (Cutler & Lindblom 1987). 
Molecular viscosity alone is likely to be less effective in 
bringing the star into uniform rotation than magnetic 
braking. Nascent neutron stars may also be subject to 
convective instabilities (e.g. Pons et at. 1999), but the role 
of convection in rotating magnetic stars is not well under- 
stood (cf. Tassoul 1978). 

For weak magnetic fields and high values of (3, the neu- 
tron star merger remnant is likely to develop a bar. The ac- 
companying quasi-periodic gravitational wave signal may 



be observable by the new generation of gravitational wave 
laser interferometers under construction (Lai & Shapiro 
1995; Shibata, Baumgarte, & Shapiro 1999b). 

For strong magnetic fields and small values of (3, mag- 
netic braking is likely to dominate the evolution of differ- 
entially rotating neutron stars, and may alter the velocity 
profile within minutes. On this timescale, differential rota- 
tion will no longer be able to support hypermassive stars 
formed in binary merger. In the resulting delayed col- 
lapse, a brief secondary burst of gravitational waves will 
be emitted. The frequency of this secondary burst may 
be quite high^, but since the angular momentum param- 
eter J/A'P may be close to unity, the amplitude could be 
large enough to be observable by an advanced generation 
of gravitational wave detectors (Stark & Piran, 1985). If 
the orbital parameters, including the masses and radii of 
the stars, can be determined during the inspiral and early 
merger phase, and if the time of the initial coalescence can 
be inferred from the initial burst signal, then the measure- 
ment of this delay in the final collapse may provide an 
estimate for the strength of the wound-up magnetic field 
in the interior of the merged neutron star. 

4. DISCUSSION 

We find that the maximum mass of a differentially ro- 
tating star can be significantly higher than that of nonro- 
tating or uniformly rotating stars, even for modest degrees 
of differential rotation. As an immediate consequence, it 



is possible that binary neutron star coalescence does not 
lead to a prompt black hole formation, but that, instead, a 
differentially rotating, hypermassive quasiequilibrium neu- 
tron star is formed. This has important consequences for 
the gravitational wave signal from such an event and pos- 
sibly for the prospects of explaining gamma ray bursts by 
binary neutron star mergers. 

Pulsars are likely to be uniformly rotating, since mag- 
netic braking and viscosity will bring any initially differ- 
entially rotating stars into uniform rotation. The well es- 
tablished maximum masses of uniformly rotating neutron 
stars are therefore relevant for old neutron stars, includ- 
ing millisecond pulsars, while our much higher maximum 
masses may be relevant for nascent neutron stars in a 
transient phase in a supernova, following fallback, or in 
a merged binary. 
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which yields / -^ 4 kHz for M = 3Mq ; 



